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ABSTRACT 

The location of warm dust producing the Mid-infrared (MIR) emission in Type 1 Active Galactic 
Nuclei (AGNs) is complex and not yet fully known. We explore this problem by studying how the MIR 
covering factor (CFmir = -^Mm/^boi) correlates with the fundamental parameters of AGN accretion 
process (such as Lboi, black hole mass -Mbh, and Eddington ratio L/L-^dd) and the properties of 
narrow emission lines (as represented by [Olll] A5007), using large data sets derived from the Sloan 
Digital Sky Spectroscopic Survey (SDSS) and the Wide Infrared Sky Survey (WISE). Firstly we find 
that the luminosity of the [O III] wing component (£ w i n g) correlates more tightly with the continuum 
luminosity (ALa(5100)) than the luminosity of the line core component (L C ore) does, which is in line 
with our previous conclusion that the wing component, generally blueshifted, originates from the polar 
outflows in the inner narrow-line region (NLR). We then find that the MIR CF shows the strongest 
correlation with £ W mg/£boi rather than with L core /Lboi or the above fundamental AGN parameters, 
and the correlation becomes stronger as the infrared wavelength increases. We also confirm the anti- 
correlations of CFmir with Lboi and Mbh, and the lack of dependence of CFmir on the Eddington 
ratio. These results suggest that a large fraction of the warm dust producing MIR emission in AGNs 
is likely embedded in polar outflows in the NLR instead of in the torus. 

Subject headings: galaxies: Seyfert — galaxies: kinematics and dynamics — infrared: galaxies — 
quasars: emission lines 



1. INTRODUCTION 

The mid-infrared (MIR) emission of radio quiet active 
galactic nuclei (AGNs) may come from the dusty torus 
and/or narrow line region (NLR), as well as star forma- 
tion regions in the host galaxies. High-resolution imag- 
ing in the near-infrared (NIR) provides direct evidence 
for the torus emission in a few nearby AGNs (Bock et 
al. 2000; Radomski et al. 2003). At longer wavelengths, 
observations showed extended MIR emission regions, co- 
incident with the inner NLR, in addition to the unre- 
solved core, presumably the optically thick torus (Tresch- 
Fienbcrg et al. 1987; Braatz et al. 1993; Cameron et al. 
1993; Packham et al. 2005; Honig et al. 2010). Besides, 
MIR spectroscopic analysis revealed an emission compo- 
nent, dominated by PAH emission, from star formation 
regions, and its contribution increases with increasing 
wavelength and decreasing AGN luminosity (e.g. Deo et 
al. 2009). 

Recent high-resolution observations with the MIDI 
instrument on VLT demonstrated that in NGC1068, 
NGC 424, the Circinus galaxy and probably also in NGC 
3783, most of the MIR emission lies along the polar di- 
rection, orthogonal to the major-axis of the NIR emission 
region (Raban et al. 2009; Honig et al. 2010, 2012; Tris- 
tram et al. 2012). Honig et al. (2012) proposed that the 
polar dust is entrained in radiatively driven dusty winds, 
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originated from the outerpart of the accretion disk (e.g. 
Konigl & Kartje 1994; Keating et al. 2012), or evapo- 
rated from the dusty torus (Dorodnitsyn et al. 2008). 
Thus the strength of MIR emission may be closely con- 
nected to the properties of the gas outflow. The outflow 
in the NLR is manifested by the blue-shifted and rel- 
atively broad component in the high-ionization narrow 
emission lines (e.g. Zhang et al. 2011). The NLR in 
nearby galaxies has been successfully modeled by biconi- 
cal outflows on scales of a few to tens of parsecs from the 
nuclei in nearby Seyfert galaxies (Heckman et al. 1981; 
Das et al. 2005; Crenshaw et al. 1999, 2010). It is well 
known that the NLR contains dust. It would be interest- 
ing to explore if the MIR properties are linked to those 
of narrow emission lines. 

The Wide-field Infrared Survey Explorer (WISE; 
Wright et al. 2010) surveyed the whole sky at 3.4, 4.6, 
12 and 22 /jm at sensitivities about 100-1000 better than 
previous all sky surveys. Taking the advantage of the 
WISE and the Sloan Digital Sky Survey (SDSS) data re- 
leases, it is possible to construct a large sample of AGNs 
with both measurements in MIR and optical to explore 
the connection between the MIR properties and the AGN 
fundamental parameters or NLR properties. In what fol- 
lows, we adopt a ACDM cosmology with Hq = 70 km 
s^Mpc" 1 , fi m = 0.3, and tt A = 0.7. 

2. THE SAMPLE AND THE MEASUREMENT 

2.1. The Sample and Spectral Fit 

From the spectroscopic data set of the SDSS Data Re- 
lease 4 (Adelman-McCarthy et al. 2006), we select 4178 
Seyfert 1 galaxies and quasars (namely type 1 AGNs) 
that have reliable measurements of emission line fluxes 
and the continuum suffers little host galaxy contamina- 
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tion. The sample is described in detail in Dong ct al. 
(2011). As described in Dong ct al. (2008), wc fit simul- 
taneously the continuum and emission lines in the spec- 
tral regimes near H/3+[OlIl] and Ha+[Sll] with a model 
consisting of several components: a power-law represent- 
ing the featureless AGN continuum, the Fell templates, 
and a number of gaussians for emission lines. Each of the 
narrow emission lines is fitted with a single Gaussian, ex- 
cept for the [Olll] AA4959, 5007 doublet lines, which are 
modeled with two Gaussians. One accounts for the line 
core and the other for a possible blue wing. The details 
of the fitting algorithm and spectral parameters can be 
found in Dong et al. (2011) and Zhang ct al. (2011). 

2.2. Mid-Infared Data 

We cross-correlate our SDSS Type 1 AGN catalog with 
the WISE All Sky Data Release (Wright et al. 2010) cat- 
alog using a match radius of 3 arcseconds to construct 
the mid-infrared detected type 1 AGN sample. Nearly all 
sources in our sample are detected by WISE (4174 out 
of 4178). We refer to the WISE Vega magnitude as Wl, 
W2, Wi and WA for the 3.4, 4.6, 12, and 22/im bands, 
respectively. We reject the sources that have signal to 
noise ratio less than 3 (S/N< 3) in any WISE band. Ra- 
dio loud AGNs are rejected according to the catalog by 
Lu et al. (2010) to avoid contamination to [O III] flux and 
continuum flux. The final sample has 3336 sources, and 
we call it 'extended sample' hereafter. We further make a 
redshift cut of z<0.3 to minimize potential cosmological 
evolution and uncertainty in the fc-correction that might 
complicate the correlation analysis. This 'base sample' 
consists of 608 sources. The magnitude in each band is 
converted to vF v according to the published WISE zero 
points, and the fc-correction is made by interpolating the 
Spectral Energy Distribution (SED) at rest-frame wave- 
lengths of 3.4, 4.6, 12, and 17/im. The fc-corrected lumi- 
nosities in each band are marked as LI, L2, L3 and LA. 
For every source we connect the adjacent bands on the 
\og\F\ versus log A diagram to approximate the MIR 
SED, and integrate the SED from 3.4 to 17/xm to obtain 
the total MIR luminosity, denoted as ^wise- 

3. DISTINCTION BETWEEN THE [O III] CORE AND WING 
COMPONENT 

Zhang et al. (2011) found that the equivalent width 
(EW) of the core component of [Olll] correlates nega- 
tively with ALa(5100) (Baldwin effect) while the wing 
component displays no or a weak anti-Baldwin effect. 
The different behaviors suggest that the two components 
come from distinct gas components regulated by differ- 
ent physical processes. The gaseous kinematics of the 
[O III] core traces the motion of stars in the gravitational 
potential well of bulges (e.g. Nelson et al. 2000; Wang & 
Lu 2001; Komossa et al. 2008), indicating that the [O III] 
core is emitted by the photo-ionized ISM in the galaxy 
bulge. This is supported by imaging spectroscopic ob- 
servations of nearby Seyfert galaxies, which show that 
the gaseous and stellar kinematics are coupled (e.g. Du- 
mas et al. 2008). On the other hand, the wing is much 
broader than the core but narrower than broad lines, in- 
dicating that its emitting region is much closer to the 
nuclei than the core component. Furthermore, the wings 
are usually blucshiftcd up to 1000 km s _1 , suggesting 
they are from an outflow component. High spatial resolu- 



tion spectroscopic observations of nearby Seyfert galaxies 
revealed individual high velocity knots up to a velocity 
of 3000 km s _1 on scales from parsecs to up a hundred 
parsecs (e.g. Crenshaw ct al 2010). All these observa- 
tions suggest that the wing component is formed in the 
outflow at a few to a few hundred parsecs from the nuclei. 

To examine the difference between the core and the 
wing in detail, we plot the luminosities of the [O III] core 
and wing against ALa(5100) in Figure 1 using the ex- 
tended sample. Obviously, both the core and the wing 
luminosities correlate with ALa(5100). The best fit of 
a log-linear model yields a slope of 0.86±0.01 for the 
core component, and 1.14±0.01 for the wing component. 
These slopes are consistent with the analysis of Zhang ct 
al (2011) on the Baldwin effect. However, the plot reveals 
an interesting feature that the correlation is much tighter 
for the wing than for the core. Qualitatively, the RMS 
scatters to the best fit of the log-linear model are 0.29 dcx 
for the wing, and 0.42 dex for the core. The difference is 
significant according to F— test with a chance probabil- 
ity P < 10~ 26 . The i W ing correlates more tightly with 
the the broad Ha luminosity, which used for bolomet- 
ric luminosity estimation in next section, than the L C oro 
does for the basic sample. Note that the decomposition 
into the wing and the core, and other measurements in- 
troduce an uncertainty of about 0.15 dex in the EWs of 
the wing and the corco 

The tighter correlation between the wing and the con- 
tinuum luminosity indicates that outflows do not change 
dramatically from object to object, while the larger scat- 
ter in the core component suggests that the variations of 
ISM states like density, density profile, ionization state 
etc., are pronounced. Physically, the ISM gas in the 
bulge is regulated by several physical processes, such as 
gas accretion, star-formation, stellar wind and AGN feed- 
back. The competition between different processes may 
lead to very diverse physical conditions of the gas among 
AGNs, which introduces a scatter in the correlation. In- 
terestingly, Zhang ct al. (2011) found that the EW of the 
[O III] core correlates strongest with the outflow velocity 
of the whole [Olll] profile (the outflow velocities being 
negative due to the outflow, hence the correlation), sug- 
gesting the interplay of the outflow and ISM of the galaxy 
bulge. This can be interpreted most directly in terms of 
AGN feedback: high velocity outflows blow away the ISM 
in the host bulge. Also, continuum variations may intro- 
duce a larger scatter in the core component than in the 
wing component because the strength of the core com- 
ponent depends on the long-term variability while the 
strength of the wing component follows more closely the 
continuum variations (Zhang ct al. 2013). The effect 
of the variability depends on the light curves over time 
scales of 10 3 years. 

4. CORRELATION BETWEEN C*F M ir AND AGN 
PROPERTIES 

We use the ratios of mid-infrared to bolometric lumi- 
nosities to quantify the MIR CFs, i.e. the relative con- 

4 For every source, we use best-fit 2-Gaussian model of the 
[O III] A5007 line plus the measurement errors of the SDSS spec- 
trum to generate a new mock spectrum, then fit the mock [OIII] 
line with two Gaussians. Statistically in the extended sample, the 
standard deviations of the difference between the output and in- 
put fluxes are 0.12 and 0.15 dex for the core and wing components, 
respectively. 
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tribution of MIR to bolomctric luminosity. The bolo- 
metric luminosity is estimated from the broad line Ha 
luminosity using a bolomctric correction factor of 170, 
which is combining Lboi=12.17x ALa(5100) in Richards 
et al. (2006) and AL A (5100)=14xL H a in Stern & Laor 
(2012)0 This is similar to the 130* £4 adopted in Stern 
& Laor (2012). We do not use ALa(5100) directly as an 
estimate for the bolomctric luminosity, because it may be 
subject to host contamination and the optical-UV slope 
gets bluer in higher luminosity AGN (See Section 4.1 in 
Zhang et al. (2013) for details). The median bolometric 
luminosity of our sample is 10 erg s" 1 . In Figure 2 we 
plot the histograms of mid-infrared luminosity to bolo- 
metric luminosity for the individual wavelengths (Panels 
a to d) and for the total MIR luminosity (Panel e). All 
histograms show a log- normal distribution. For CFwise, 
the distribution is centered at log CFwise = -0.37 (43%) 
and has a dispersion of 0.243 dex. To first order, this 
means a covering factor of 43%, and is consistent with 
the type 2 fraction of -~50% at this luminosity in liter- 
ature (e.g. Hao et al. 2005; Lawrence & Elvis 2010). 
However, one needs to caution that the type 2 fraction 
may not be equal to the covering factor of warm dust. 
Anisotropy in primary continuum and IR emission, pres- 
ence of large scale obscuring material and hot dust that 
emits at wavelength shorter than 3.4pm lead to deviation 
of one from another (Calderone et al. 2012; Roseboom 
et al. 2013; Ma & Wang 2013). 

The CFmir spreads over a dynamic range (the range 
it spans) of more than 1.8 dex at the 17pm band and 
the dynamic range increases with increasing wavelength. 
This trend is clearly seen in Figure 2 in panel (a) to (d) 
respectively. The cr's of CFmir increase from 0.200 dex 
at 3.4/im to 0.298 dex at 17pm. One possible reason 
of increasing dispersion with wavelength is that star- 
formation contributes a large fraction of emission at 
longer wavelengths. As discussed at the beginning of Sec- 
tion 5, star-formation would contribute less than 10% at 
17/im, thus introducing a dispersion less than 0.04 dex. 
This alone is not enough to account for the 0.22 dex 
difference between dispersions of CFl and CF4. The 
intrinsic variance after subtracting the measurement er- 
rors in IR luminosity still show a difference of 0.20 dex. 
There are two possibilities for this difference: 1) The 
17 pm emission is from the dust torus but the fraction 
of warm dust in the torus varies drastically from object 
to object due to either different degree of dumpiness or 
different optical depth of the torus. 2) There is an addi- 
tional warm dust component, and the contribution from 
the additional component varies. As discussed in Sec- 
tion 5.1, the torus explanation is inconsistent with our 
findings, and we argue below that the latter is plausible. 

4.1. Correlations between CFyim and fundamental 
AGN properties 

First, we analyze how CFmir correlate with redshift 
and the AGN fundamental parameters, such as Lboi, 
Mbh and L/Lsdd using the basic sample. The results 
are summarized in column 1, 2, 3, 4 of Table 1. We 
calculate the black hole masses based on the FWHM 

5 The bolometric correction may depend on the accretion rate 
and luminosity (e.g. Vasudevan & Fabian 2007; Lusso et al. 2012). 



of broad H/3 and AL,\(5100) using the formalism pre- 
sented in Wang et al. (2009, their Eqn. 11). The 
Eddington ratio (L/LEdd) is the ratios between the 
bolometric and Eddington luminosities: L/Lsdd = L bQl > 
L E dd = 1-26 x 10 38 (A/bh/Mq) erg s" 1 . We use 
Spearman correlation coefficient (p) as a non-parametric 
measure of the statistical dependence between two 
quantities. 

We find that the CF M ir correlate negatively with 
Lboi and Mbh, but very weakly with L/L^dd- We get 
CFl oc L bo r 21 , CF2 cx L bo r a21 , CF3 oc L bo r°' 3 \ 
CF4 oc L bo i °- 35 , and C^wise/ oc Lboi 26 respectively 
as shown in Figure 3. The typical error of the slopes 
is 0.03. These results are fully consistent with previ- 
ous works (Maiolino et al. 2007; Treister et al. 2008; 
Calderone et al. 2012; Ma & Wang 2013; Roseboom et al. 
2013). We plot the best regression fits in Calderone et al. 
(2012) and Ma & Wang (2013) on Panel (e) to illustrate 
the consistency. Calderone et al. (2012) integrated the 
MIR luminosity at 2. 1-13. 7pm, and Ma & Wang (2013) 
used 3-10pm luminosity. To make a fair comparison, we 
add 0.184 dex to the CF Mm in Ma & Wang (2013) to 
account for the shorter wavelength range. The 0.184 dex 
is the mean ratio of 3.4-17pm luminosity to 3. 4-11. 6pm 
luminosity in our sample. We do not make any correction 
to the regression fit in Calderone et al. (2012) because 
the gain in 2-3. 4pm approximately equals to the loss in 
13.7-17pm according to the mean SED of all quasars in 
Richards et al. (2006). We can see that the fitting lines 
from literatures and our best regression fit are consistent 
with each other at this luminosity range. This further 
supports the bolomctric correction we adopt is reason- 
able. 

The slopes of CFmir versus Mbh relations are -0.22, 
-0.23, -0.35, -0.42, -0.30 for LI, L2, L3, L4, and L W ise 
respectively, consistent with previous works (e.g. Mor et 
al. 2011; Ma & Wang 2013). The lack of dependence of 
CFmir on L/L^dd is reported by many authors (Cao et 
al. 2005; Mor et al. 2011; Ma & Wang 2013). Consider- 
ing the limited range of L/iEdd(l-5 dex) for our sample, 
we refrain from over- interpreting this result. 

4.2. Correlations between CFmir and narrow emission 
line properties 

Next, we explore whether the MIR emission is linked to 
the properties of optical narrow emission lines as repre- 
sented by [O III] A5007. We present the results of correla- 
tion analysis between CFmir and the [O III] total, core, 
wing luminosity to Lboi ratios using the basic sample. 
The results are summarized in column 5, 6, 7 of Table 1. 
The CFmir in all 4 bands correlate significantly with 
^wing/iboi- Among them, the strongest one is between 
CF4 and L w i ng /L bo i with p = 0.52. The correlation 
weakens as wavelength decreases: with p =0.47, 0.36, 
0.32 between L wing /L bo i and CF3, CF2 and CFl, re- 
spectively. This is in parallel with the increasing dynamic 
range of the CFs with wavelengths. Note that the bolo- 
metric correction has RMS scatters of order 0.38 dex. 
So it is likely that the weakening of the correlations is a 
combination of decreasing dynamic range and relatively 
large uncertainty in the bolomctric correction. The phys- 
ical reason of increasing dynamic range with increasing 
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wavelength is discussed in Section 5.3. We plot CFA ver- 
sus I/wing/^boi in Figure 4. A log linear fit yields a slope 
of 1.19±0.04 with an RMS scatter of 0.308 dex. This 
scatter is significantly larger than the combinations of 
measurement errors: 0.11 dex in L4/Ljja and 0.15 dex 
in L w i ns /Ln a , indicating that there is a considerable in- 
trinsic scatter in this relation. No correlation is found 
between outflow velocity and any IR properties in our 
sample. 

5. DISCUSSION 

To facilitate discussions on the implication of the cor- 
relation between L w i n g/^boi and L4/Lboi, we note that 
the star-formation contamination is likely small in the 
MIR band for these luminous type 1 AGNs in our sam- 
ple. For most of our sources, the infrared luminosity at 
17/xm is in the range of 10 43 to 10 46 erg s _1 (10 9 ' 5 to 
10 12 ' 5 L©). According to Fig 10 of Deo et al. (2009), 
starburst contribution is likely less than 10% at 17/xm 
and a few percent in the LI band. This is more likely to 
be true because we have removed any source with sub- 
stantial stellar contamination in optical. In the following 
sections, we will discuss the location of the warm dust 
and its relation with torus in light of the above results. 

5.1. Classical Steady Torus? 

Judging from studies of AGN at subarcsecond resolu- 
tion in the MIR, the majority of MIR emission originates 
from the central few pc region (Jaffc et al. 2004; Tris- 
tram et al. 2007). In NGC 1068, about 75% of the MIR 
emission is from the central unresolved core (Bock et al. 
2000). This was believed to be the torus, and a variety of 
torus models were constructed to match the NIR-MIR 
SED (e.g. Fritz et al. 2006; Honig et al. 2006; Schart- 
mann ct al. 2008; Nenkova et al. 2008a, b). If the dust 
dominates the MIR emission, as the torus acts to confine 
the ionization cone, there should be an anti-correlation 
between the torus subtending angle and the covering fac- 
tor of NLR; we thus should obtain a negative correlation 
between CF4 and L COIC / 'Lboi- This is however in contra- 
diction to the observation. 

5.2. Dusty NLR clouds? 

NLR is theoretically expected to contain dust (Net- 
zer & Laor 1993; Dopita ct al. 2002). Obscrvationally, 
extended polar MIR emission is found around the com- 
pact cores of several AGNs (See Section 1). Groves et al. 
(2006) modeled the narrow line spectrum of NGC 4151 
with dusty photo-ionization models, and found that 10- 
50% of its 25/xm emission is from the NLR dust. The 
imaging observations also indicate that >27% of MIR 
emission emerges from a region coincident with the NLR 
in NGC 4151 (Radomski et al. 2003). The warm dust ex- 
tends from the unresolved core to about lOOpc. Mor et al. 
(2009) fitted the MIR spectra of PG QSOs with 3 compo- 
nents: a black body, the clumpy torus emission (Nenkova 
et al. 2008a,b) and the NLR dust. They derived that 
above 18/xm, the NLR dust contributes as much as 40% 
of the emission. Schweitzer et al. (2008) demonstrated 
the the NLR dust could produce the silicate emission 
features observed in AGNs. As the cold gas and the 
dust are coupled in the NLR, this would naturally ex- 
plain the correlation between CF4 and i w i ng /iboi and 



the correlation weakens at shorter infrared wavebands, 
qualitatively. 

However, the NLR dust models do not consider the 
dust in the outflow component of NLR, so they can not 
explain why CFA is much more strongly correlated with 
-^wing/iboi than with L core /Lboi- Furthermore, most of 
the infrared emission from 5 to 18/xm is attributed to 
the torus emission, which is not consistent with more 
recent higher resolution imaging observations in several 
nearby Seyfert galaxies (see §1). In Mor ct al. (2009), 
the distance of NLR dust is estimated to be Rnlr = 
295 x L®q pc, where L4Q is the bolometric luminosity in 
unit of 10 46 erg s _1 . In high luminosity sources, we lack 
observational evidence that the outflow may reach that 
far from the nuclei. Although detailed model is needed to 
predict the infrared spectrum, it is clear that any optical 
thin dust in the outflow region will be hotter than that 
in the more distant classical NLR, and will contribute 
significantly to emission below 15 /xm. Alternatively, the 
dust clouds in outflows are optically thick to infrared 
light, radiative transfer effect would result in much colder 
infrared emission. We therefore do not think that dusty 
NLR clouds on scales of >100 pc are responsible for our 
observed correlation of L w - lng / Xboi with CFmir- 

5.3. Dusty Outflows? 

It was pointed out by Krolik and Begelman (1988) that 
a steady doughnut torus suffers the problem that there 
is not enough vertical motions capable of sustaining the 
H/R~l geometry. One solution to this problem is assum- 
ing that the dust is entrained in the magnetically driven 
accretion disk wind (Emmcring, Blandford & Shlosman 
1992; Elitzur & Shlosman 2006), and the outer part of 
the outflow acts as a clumpy torus. Alternatively, the IR 
radiation pressure on dust may also offer enough pres- 
sure to support a thick torus, even launches an outflow 
(Krolik 2007; Dorodnitsyn et al. 2011, 2012). In this 
type of dynamical models, the outflow and the torus arc 
naturally linked. A connection between the outflow and 
the dust emission is expected. The detection of dominant 
MIR emission in polar direction led Honig et al. (2012) to 
propose that the polar dust is entrained in a radiatively- 
driven dusty winds. This dust component is an inward 
extent of NLR dust observed in many nearby AGNs but 
lies only a few pc away from the nuclei. This is where 
the outflow originates and accelerated. In this paradigm, 
the outflow gas and the warm dust share similar cover- 
ing factor, and would lead to our finding of a positive 
correlation between CF4 and £ w j ng /Lboi naturally. The 
contribution of the polar warm dust to the mid-infrared 
emission increases with wavelength, so the correlation 
between C-Fmir and L w - lng / Lboi strengthens with wave- 
length. By comparison, the static classical NLR traced 
by the [O III] core component lies further from the nuclei. 
The dust embedded in the static NLR emit dominantly 
at wavelength longer than 20^m, thus only connects to 
the warm dust loosely. If the correlation between the 
warm dust and the outflow strength is really driven by 
the dust in the polar winds, an implication is that the pc 
scale MIR elongation in NGC 1068, the Circinus galaxy 
and NGC 424 may not be special cases, but a common 
phenomenon for Seyfert galaxies. Most of the warm dust 
may reside in the outflow instead of in the torus. Besides, 
the outflow may have a significant effect on the dust re- 
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distribution by transporting the dust from the central 
region to larger scales. 
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TABLE 1 

Correlation Tables between the IR covering factors and the optical parameters a 





z 




£bol 




M BH 




L/ L E dd 


L[o in] / Lboi 


icoro/Lbol 


-^wing/-^bol 


C-Fwise 


0.06(lc 


■01) 


-0.40(3c 


25) 


-0.38(2c 


22) 


-0.10(7c-03) 


0. 28(lo-12) 


0.09(2c-02) 


0.45(lo-32) 


CFl 


0.08(2e 


02) 


-0.37(5e 


22) 


-0.32(8e 


17) 


-0.12(lo-03) 


0.10(le-02) 


-0.07(6e-02) 


0.32(3e-16) 


CF2 


0.09(le 


1)2) 


-0.35(4e 


20) 


-0.33(2e 


17) 


-0.09(le-02) 


0.14(5e-04) 


-0.04(2c-01) 


0.36(8e-20) 


CF3 


0.05(le 


01) 


-0.40(4e 


25) 


-0.38(5e 


23 ) 


-0.10(lc-02) 


0.31(2c-15) 


0.12(le-03) 


0.47(4e-35) 


CF4 


0.02(5e 


01) 


-0.42(2c 


27) 


-0.40(3c 


25) 


-0.10(6c-03) 


0.43(9e-29) 


0.25(9e-ll) 


0.52(lc-43) 



a For each entry, we list the Spearman rank correlation coefficient (p) and the probability of the null hypothesis that the 
correlation is not present (P n uii) in parenthesis, for the 608 objects in the basic sample. p >0.50 are marked in Boldface. 
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Fig. 1. — Left: The [OIII] core luminosity versus AL^(5100) using the extended sample (3336 sources). The solid line shows the best 
fitting line, and the cyan contours show the density of the distribution. The levels of the contours are 15, 30, 45, 60 and 75 respectively. 
Right: The [OIII] wing luminosity against AL^(5100) using the extended sample. The levels of the contours are 15, 30, 60 and 80. 
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Fig. 2. — Distributions of MIR covering factors: C-Fmir = LMIR/Lbol using the basic sample. The solid line is the gaussian fit to the 
distribution, the centroid and a of the gaussian are shown in the left-up corner. Panel (a) to (d) are C-Fmir distributions at each band, 
and Panel (e) is the distribution of CFwise = ^wiSE/^bol- 
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Fig. 3. — C-Fmir against Fbol °f the basic sample. The blue line is the best regression fit, and the red line is the mean value in each Lboi 
bin with binsize of 0.3 dex starting from Lboi = 10 43 ' 8 erg s —1 . The slope of the regression fit is shown in the left-down corner. Panel (a) 
to (d) are CFmir versus Lboi at each band, and Panel (e) is the CFwise versus the bolometric luminosity. The green and purple lines are 
the results from Ma & Wang (2013) and Calderone et al. (2012). 




Fig. 4.— OF4 versus _L w i n g/Z/k i using the basic sample (608 sources). The correlation coefficient is shown in the left-up corner. The 
solid line is the best regression fit. 



